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Summary 

Equivalent crystal theory (ECT) is a new semiempirical 
approach to calculating the energetics of solids with defects. 
The theory has successfully reproduced surface energies in 
metals and semiconductors. The theory of binary alloys, both 
with first-principle and semiempirical models, to date has not 
been extremely successful in predicting the energetics of 
alloys. This paper presents an extension of ECT, which is used 
to predict the heats of formation, cohesive energies, and lattice 
parameters of binary alloys of Cu, Ni, Al, Ag, Au, Pd, and 
Pt as functions of composition. This procedure accurately 
reproduces the heats of formation versus composition curves 
for a variety of binary alloys. The results are then compared 
with other approaches, such as the embedded atom method 
or Miedema’s method. In addition, a new sum rule, which 
predicts the cohesive energies and lattice parameters of alloys 
from pure metal properties more accurately than Vegard's law, 
is presented. 


Introduction 

The equivalent crystal theory (ECT) (refs. I to 3), recently 
extended to include a proper treatment of bond length 
anisotropy and bond angle anisotropy (ref. 4), is still limited 
in its inability to treat systems with more than one atomic 
species. The basic idea underlying this theory, that is, the 
existence of an equivalent, fictitious crystal for each atom 
encompassing the range of a certain local defect, is not appli- 
cable to alloys, since the only difference between this fictitious 
single crystal and the actual ground state crystal is that its lattice 
constant differs from its ground state value. Thus, nothing in 
ECT accounts for a variety of atomic species in the local 
environment. In other words, an alloy cannot be considered 
as a defect of an otherwise pure crystal (i.e., a crystal of one 
atomic species). 

Retaining the atomic identity of the constituents is essential 
in developing a technique for calculating alloy properties. 
Without atomic identity, detailed knowledge of the particular 
alloy properties we wish to predict would be required. Such 
particularity would limit the range of applicability of any 
method attempting to deal with multicomponent systems. To 
use the current version of ECT for alloys, it would be 
necessary to redefine the system, its components, and the 
formal perturbative series so that each atom would become 


an “average” alloy particle and so that the interactions with 
neighboring alloy particles would be redefined accordingly. 
However, to apply ECT directly, a priori knowledge is 
required of the structure and properties of that alloy with the 
lowest possible cohesive energy, which is clearly impossible 
as in most cases it is one of the pure components, and not the 
alloy, that is sufficiently well known. The values of the binding 
energy and cohesive energies are the same at the minimum 
in the binding energy curve when it is referenced to zero at 
infinite separation. 

In deriving an alternative approach to an exact treatment 
of alloys (that, if satisfied, would allow for a simple and 
efficient study of alloys and defects), we face two constraints: 
( 1 ) We must keep a single-species description, thus allowing 
us to treat individual atoms as building blocks of any metallic 
compound, and (2) we must be able to introduce specific alloy 
properties in a perturbative fashion to properly account for 
the behavior of multicomponent systems. Further, in order 
to retain the simplicity and numerical accuracy of the original 
ECT, we chose to base our approach on it. 

In the original ECT, each atom in a defect crystal is assigned 
an equivalent crystal of the same atomic species. The lattice 
parameter of this equivalent crystal is determined such that 
the energy of an atom in the equivalent crystal is the same 
as the energy of the atom in the defect crystal. The lattice 
parameter of the equivalent crystal is obtained via a per- 
turbative scheme which translates into solving a simple 
transcendental equation containing information about the 
nature of the defect. Once the equivalent lattice parameter is 
determined, the energy is computed by means of the universal 
binding energy relation. This energy is actually the difference 
in the energy of the atom in the defect crystal and that in the 
atom in the ground state crystal. 

In an alloy, however, an atom of a given species finds itself 
in a different environment from the one in the ground-state 
pure crystal of its own species: the geometry is different, and 
some of the neighboring atoms are of a different species. These 
two changes should be dealt with separately, as they cannot 
both be considered defects. Within the framework of ECT, 
a defect is any change of the local environment of the pure 
crystal where nothing is necessarily conserved, except the 
identity of the atoms. The equivalent crystal of any given atom 
in the defect crystal is a compressed or expanded version of 
the ground-state pure crystal. In this sense, then, the formation 
of an alloy A-B cannot be considered as a “defect A” or 
“defect B” crystal. 
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The basic difference between the changes in geometry 
related to alloy formation (with reference to pure single 
crystals) and the composition effects, forces us to establish 
a procedure that will treat each aspect separately. 

We now formulate a technique that maintains the simplicity 
of ECT and uses pure metal properties to calculate alloy 
properties. This procedure is then tested by comparing its 
predictions with experimental data, first-principles calcula- 
tions, and other semiempirical and empirically based methods. 

In the section on equivalent crystal theory we briefly review 
the fundamental ideas of the equivalent crystal theory for pure 
metals. In the section 4 ‘Formalism' ’ a detailed description of 
the new formalism is given. The following two sections give 
the results and an extended discussion which includes 
comparisons with other methods. 
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equilibrium lattice constant, A 
equilibrium lattice constant 
equilibrium bulk modulus, GPa 

ratio between nearest (next-nearest) neighbor 
distance and the lattice parameter 

energy 

heats of solution of impurity in host 
cohesive energy, eV 

minimum total energy, i.e., the equilibrium 
binding energy, eV 

excess internal energy for disordered alloys 
binding energy of ordered alloy 
excess internal energy for ordered alloys 
binding energy curve for (dis)ordered compound 
chemical energy (eq. (9)) 
strain energy 

heat of formation, the difference in energy 
between the constituent atoms of the compound 
at equilibrium and their composition-averaged 
bulk values in the pure crystals 

identity of the neighbor of atom / located at a 
distance r, from the host atom 

screening length, = (E c / 127rZ?r wse ) ,/2 , A 
scaling length 

nearest (next-nearest) neighbors 

number of nearest and next nearest neighbors in 
a perfect (equivalent) crystal 


n principal atomic quantum number 

P(x) interpolating polynomial for heat of formation 

of disordered alloys 


p =2/7-2 


Q 


R\ (R 2 ) 


r 

n 

^*wse 

X 
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&AB’ &BA 

A 00 

L -*m 

x 

v„(r) 
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arbitrary point on the binding energy curve of 
a certain crystal 

nearest (next-nearest) neighbor distance in the 
perfect crystal (eq. (13)) 

lattice parameter 

distance between y lh neighbor and atom under 
consideration 

equilibrium value of Wigner-Seitz radius 
arbitrary concentration 

parameter that will primarily reflect the structure 
of the electron density in the overlap region, 
parameter computed by requiring agreement 
with the experimental vacancy formation 
energy 

perturbative parameters 
energy reference 
electronic screening length 
many -body interaction potentials 

multisite correlation functions defined on an 
/? th order cluster 


Subscripts: 

A 

AB 

B 

BA 

m 

v 

x 


metal A 

impurity A in host B 
metal B 

impurity B in host A 
indicates ordered alloys 
virtual crystal 
indicates disordered alloys 


Superscripts: 

A metal A 

AB impurity A in host B 

B metal B 

BA impurity B in host A 


Equivalent Crystal Theory of Metal 
and Semiconductors 

The equivalent crystal theory (ref. 4) is a new tool for the 
treatment of real material defects at the atomic level. The 
method treats surface energies and the surface relaxation of 
semiconductors and metals accurately. The basic idea of this 
method is that, for any crystal with a defect, an atom in the 
vicinity of the defect has the same energy that it would have 
in a certain perfect (equivalent) crystal, which we define later. 
The defect crystal can be formally described as a perturbation 
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Scaled lattice constant, a* 

Figure 1.— Scaled binding energy per atom of a crystal as function of inter- 
atomic separation for representative solids. The sources for unsealed results 
are listed in ref. 12. 



Scaled lattice constant, a* 

Figure 2.— Scaled binding energy versus scaled separation for representative 
cases of cohesion, bimetallic adhesion, chemisorption, and diatomic 
molecule. Sources of unsealed results are listed in ref. 12. 


of a perfect crystal whose lattice parameter is chosen to 
minimize the perturbation. 

The procedure is very simple, and for each atom near a 
defect, it involves solving a few transcendental equations that 
represent the perturbation between the real defect crystal and 
the perturbed, equivalent crystal. The number of equations 
depends on the complexity of the defect. The solution of these 
equations gives the lattice parameter of the equivalent crystal. 
In a straightforward calculation, each atom is then assigned 
its own equivalent crystal whose energetics follow a universal 
behavior given by the universal binding energy relation. 

Total energies as a function of interatomic spacings have 
been discovered to have a single, universal form for bimetallic 
adhesion, for cohesion in metals (fig. 1), for metallic and 
covalent bonds in chemisorption, for many diatomic molecules 
(fig. 2), and even for nuclear matter. All these curves are 
scaled on to one universal form, which can be obtained by 
a simple scaling of the total energy: 

E(a) = E c E*(a) (1) 

where 


and where E ( is the minimum value of the total energy, that 
is, the equilibrium binding energy, and / is a conveniently 
defined scaling length. As shown in figures 1 and 2, a simple 


analytic form accurately represents the universal energy 
relation: 

E*(a‘) = - (1 + a*) exp - a’ (3) 

From the point of view of an individual atom in a certain 
crystal, each point of the binding energy curve can be 
interpreted in two different ways. To fix ideas, let’s consider 
an arbitrary point Q on the binding energy curve of a certain 
crystal (fig. 3): 



Figure 3.— Any point Q on binding energy curve. 
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Point Q then denotes (1) the energy Eq of the atom when the 
perfect crystal is uniformly expanded until its lattice constant 
has the value gq, and (2) the energy Eq of the atom in a 
defect crystal, but where the defect is such that its equivalent 
crystal has the same lattice constant ciq. In the process of 
creating the defect, that particular atom increased its energy 
by an amount 


hE = E c - E q (4) 

The sum of similar contributions from the other atoms 
surrounding the defect amounts to the total energy necessary 
to create the defect. 

In order to determine the equivalent lattice constant ciq , we 
solve a simple transcendental equation where the only input 
is given by the atomic positions of the surrounding neighbors, 
thus carrying all the relevant information about the nature of 
the defect. This equation, which in what follows will be 
referred to as the ECT perturbation equation, results from a 
simple parameterization of an exact perturbative treatment of 
the problem: 


N\(c\a Q ) p exp (-ac\,a Q ) + N 2 (c 2 a Q ) p 



where r. stands for the distance between the j lh neighbor and 
the atom under consideration; N\ and N 2 are the number of 
nearest and next-nearest neighbors in a perfect (equivalent) 
crystal: p = In — 2, where n is the principal atomic quantum 
number; X is the electronic screening length; a is a parameter 
computed by requiring agreement with the experimental vacancy 
formation energy; and C\ (c 2 ) is the ratio between the nearest 
(next-nearest) neighbor distance and the lattice parameter. The 
sum on the right side of equation (5) is over the neighbors 
of the atom in the defect crystal. Once equation (5) is solved 
with respect to the equivalent lattice constant a Q , the change 
in energy is found from the universal binding energy relation: 

E = E ( .[l - (1 + a' Q ) exp (-a’ Q j\ (7) 


where 


a n = 


C\(a Q - a c ) 


E c is the cohesive energy, and C\ is the ratio between the 
equilibrium Wigner-Seitz radius and the lattice parameter. This 
procedure is performed for each atom that defines the defect 
and the sum of these energies represents the energy of forming 
the defect. 

A complete formulation of ECT can be found in reference 4, 
where additional terms dealing with bond length anisotropy, 
bond angle anisotropy, and face diagonal anisotropy are 
included. Including these terms does not change the essence 
of the method, which is based on the existence of an equivalent 
crystal for which the brief summary included here suffices. 
For simplicity, we will not refer to these terms when dealing 
with the extension of ECT to alloys, but it must be understood 
that in any practical application, a full ECT treatment of the 
defect should be carried out following the prescription in 
reference 4. 


Formalism 

General Considerations 

In order to fix these ideas, we will describe the ensuing 
formalism with reference to the following example. Consider 
two pure single crystals: one of atomic species A (lattice 
parameter a*) and one of atomic species B (lattice parameter 
a f). This will be the initial state. The final state will be a 
certain alloy AB with lattice parameter a x . The ideal process 
of alloy formation is shown in figure 4. 

Let us focus on one of the atoms in the A crystal. Figure 5 
represents the transformation undergone by this atom. There 
is a change in geometry (the lattice parameter changed from 
a * to a x ) and a change in composition (some of the neighbors 
are changed to B atoms, denoted by dots). As discussed above, 
we will break up this transformation into two independent 
transformations, as shown in figure 6. 

In the first transformation, the identity of the atoms is 
conserved. The atom in question (denoted in figs. 5 and 6 
by (x)) sees its environment changed only in terms of the 
relative distances of the atoms surrounding it. This is a defect 
that can be straightforwardly treated with the current ECT 



Figure 4.— Ideal process of alloy formation. 
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Figure 5.— Transformation of atom A in crystal. 
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Figure 6.— Breakup of transformation shown in figure 5 into two independent 
transformations. 


• a E - a? 

°A = — r - 

l A 

s A = + 1 if a A > 0, = — 1 if a A < 0, and a E is obtained 

by solving a transcendental equation that will be discussed 
later. 

The last term in equation (9) is introduced to compensate 
for the fact that the chemical energy, as defined here, is 
independent of the structure of the crystal. Without it, the 
chemical energy would be a constant, solely dependent on 
chemical composition, regardless of the relative positions of 
the atoms. It will be shown later that we must include this term 
in order to properly account for the asymptotic behavior of 
the enthalpies of formation of alloys. However, it will also 
be shown that for the description of equilibrium properties of 
the alloys, there is a certain arbitrariness concerning the 
definition of this term as well as its need altogether. 

Strain Energy 

The calculation of the strain energy is straightforward, as 
it only requires the knowledge of the location of the occupied 
sites surrounding the atom being studied. The first step 
involves solving the ECT perturbation equation for that atom 
in order to determine the equivalent lattice parameter a : 


without reference to alloy formation. The atom suffers a 
change in energy, e*, which we will call from now on strain 
energy, because it is related only to lattice deformations. The 
change in energy is simply 

ef = E*[ 1 — (1 + a s ) exp (-«*)] (8) 


where 


NR P \ exp (-a/?,) + Af/K exp 



= ^ r , exp [ - (a + S(r,) )r] (10) 

i 


where 


ci s = 


a — a 




COS 


7 T(r- R\’) 

(R ( 2 - R") 


( 11 ) 


and where a is the lattice parameter of the equivalent crystal 
of the atom in question, which can be easily obtained by 
applying single-crystal ECT. The parameters E* , and l A 
are the cohesive energy, equilibrium lattice parameter, and 
scaling length, respectively. 

In the second transformation, the geometry of the equilib- 
rium crystal is conserved. The atom in question sees its 
environment changed only in terms of the identity of its 
neighbors. This change in energy, which we call the chemical 
energy, is given by 


for 

R‘\ = C\a e and R 2 = c 2 a e (12) 

where /*, is the distance between the / th neighbor and the atom 
under consideration and /V (A/) is the number of (next) nearest 
neighbors in the perfect crystal. Also, 

/?, = C\a and R 2 = c 2 a (13) 


e* = s a E?[ 1 - (1 + a A ) exp (-a\)] exp ( -a A ) (9) 


where cj (c 2 ) is such that R \ (R 2 ) is the (next) nearest 

neighbor distance in the perfect crystal. Once the equivalent 
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lattice parameter a is known, the strain energy is easily 
computed using equation (8). 

Chemical Energy 

In the single-crystal ECT, where all the atoms are of the 
same atomic species, we apply perturbation theory in order 
to find the energy of the defect crystal. The perturbation is 
basically due to the difference in potentials between the defect 
solid and the ground-state crystal. As described in the original 
formulation of the method, it is reasonable to parameterize 
the first-order contributions to the perturbation expansion as 

A E oc Rf 7 exp( - aR) ( 14) 

where p = 2 n — 2 (where n is the atom principal quantum 
number) and a is a parameter that will primarily reflect the 
structure of the electron density in the overlap region. In single- 
crystal ECT the parameter a is determined for metals so that 
the energy to form a rigid (or unrelaxed) vacancy is equal to 
the experimental value. 

To a good approximation, these concepts should remain 
valid in the case of alloys, and we will adopt the same func- 
tional form used in equation (6) to describe the perturbation 
due to the dissimilar atomic species. In order to deal with 
arbitrary defects and structures in future applications, as well 
as with multicomponent systems, it is convenient to “localize" 
this effect and assume that the global property parameterized 
by a (i.e., the tails of the overlapping electron densities) can 
be separated into pairs of interacting atoms. In this approxi- 
mation, the electron density in the region between two atoms 
of the same species would be unaffected by the presence of 
neighboring atoms of a different species. The perturbation 
would then be localized in the region between two dissimilar 
atoms. This assumption enables us to define the parameter 
a AB as 

Oi ab = a a &ba ( 15 ) 


a homogeneous crystal where its constituents are quasiatoms 
that properly model the chemical composition of the impure 
crystal, then any alteration in the composition could be studied 
by recourse to a single set of equivalent crystals of the 
quasiatoms. This requires the a priori knowledge of the 
properties of this quasicrystal; however, these are the 
properties that we wish to predict. Further, the quasiatom 
approach makes the calculation of the energetics of an arbitrary 
defect, both in composition and geometry, impractical. 
Therefore, we are attempting to develop a method that retains 
the simplicity of the original ECT, where the identity of the 
individual constituents is maintained. We must, then, work 
with the framework of single crystals, where an atom of certain 
species is only “allowed" to have equivalent crystals of its 
own species, whether it is a pure environment of atoms of the 
same element or not. 

To distinguish the type of equivalent crystals one determines 
when computing strain energies, we will call these new, 
approximate equivalent crystals, which represent solely the 
changes in atomic compositions, “virtual crystals.” We must, 
of course, remember that that they are just a tool to help us 
approximate the “jump" into the otherwise inaccessible family 
of curves of the real equivalent crystals of the previously 
defined quasiatoms. 

In order to compute the net change in chemical energy 
(|? ( .|) in such a transformation (see fig. 6(b)), we first find the 
virtual crystal of the atom under consideration, where the 
defect crystal is the equilibrium crystal whose constituents are 
a mixture of atoms of different species. 

As an example, let us consider an atom of species A which, 
in its reference state, would be located in a lattice site with 
N (M) nearest (next-nearest) neighbors. The change in 
composition to that of the alloy can be represented as the 
flipping of some of the surrounding A atoms into B, C,..., 
etc., atoms. When writing the ECT perturbation equation for 
the host atom, we distinguish each interaction by means of 
the coefficients a A j(j = B,C,...), such that 


where a A is the a value for the pure metal A and A ba is a 
correction introduced by the presence of a neighbor of 
species B. The use of a AB in the ECT will be described later. 
Obviously, A ab = 0 if A = B. 

The “perturbation" parameters A.,# and A ba are the only 
new parameters introduced in this theory of alloys, and they 
will be determined by fitting to appropriate experimental data. 

The main ingredient in the calculation of e c (see eq. (9)) 
is A/., the equivalent lattice parameter for the atom embedded 
in the alloy. As noted earlier, the concept of defect, as defined 
in the ECT framework, is not applicable in this case: even 
when a single impurity is introduced in an otherwise perfect 
crystal, every single neighboring host atom has access to a 
range of atomic states that cannot necessarily be reproduced 
by an equivalent perfect crystal of a single species. If we were 
able to map the perfect crystal with the single impurity onto 


NR^ exp (— a A R^) -F MR% A exp 



= rf' exp - [a Aj + S(r,)]r, 


j 


for 


(16) 


R\ = C\ a E an d Rj — c 2 a E 

where j denotes the identity of the neighbor of atom / located 
at a distance r } of the host atom and where 

n oc af (17) 
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Once a E is known, the net change in energy is obtained by 
means of the universal expression: 

e‘‘ = s„£"[ 1 - ( I + a ' A ) exp ( - a \ ) ] exp ( - a l) (18) 

where 
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A v 
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0 X 

„A 

0 x 
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► a e 
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X X 
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ci A = 


ci t: - a e 


Figure 7.— Perfect crystal A undergoing ideal transformation. 


with s A = + 1 if a A > 0 and s A = — 1 if a A < 0, where a * is 
the equilibrium lattice parameter of a pure A crystal. 

The origin of the sign s A , which obviously accounts for the 
direction of the energy change in the alloy, can be explained 
with reference to a simple model of vacancy formation, as 
follows: Consider a perfect A crystal which undergoes an ideal 
transformation where a vacancy is created without relaxation 
of the atomic positions of the neighboring atoms. This ideal 
process is shown in figure 7, where the circle denotes the 
vacancy site and where .t* denotes one of its nearest neighbors. 

Each ® atom increases its energy by an amount e v , which 
can be easily computed with ECT. The equivalent crystal of 
each (x) atom (lattice parameter a v ), is then an expanded 
version of the ground-state crystal (lattice parameter a*). 
This process, in terms of equivalent crystals, is shown in 
figure 8. The ideal process of creating a vacancy can be 
interpreted as the uniform expansion of the corresponding 
equivalent crystals of each one of the neighboring atoms. And 
the total change in energy of the defect crystal is just the sum 
of the energies involved in each one of these expansions. 

Consider now a second, ideal, process, which is just the 
reverse of the process described before: that is, the vacancy 
is now filled, returning the defect crystal to its original state 
(as shown in fig. 9). 

Regardless of the intermediate states reached during this 
process, the net result is the lowering of the energy of the 
system to its original value. In terms of equivalent crystals, 
the equivalent crystal of each neighboring atom of the vacancy 
undergoes a uniform compression until the lattice parameter 
of the equivalent crystal exactly coincides with the equilibrium 
lattice parameter. If we focus our attention on the initial and 
final states ((1) and (3) in fig. 9), there is no net change in 
the equivalent lattice parameter, which amounts to no net 
change in the energy of the crystal. The process of filling the 
vacancy, then, is a process that lowers the energy by an amount 
proportional to the reduction in the equivalent lattice 
parameter. If it reduces to the equilibrium value, the net change 
is zero. 

As noted before, no exact statement can be made (in terms 
of equivalent crystals) with respect to a similar process where 
the vacancy is filled with an impurity atom, as shown in figure 
10. The equivalent crystal of (x) in the final state (3) can be 
taken just as an approximation to the real final state (upper 
line in fig. 10). The process of filling the vacancy with an atom 
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Figure 8. — Ideal process of tilling a vacancy. 
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Figure 10.— Ideal process where vacancy is tilled with impurity atom. 


(of the species of the host or any other atomic species) is then 
a process which lowers the energy of the original crystal. If 
the lattice parameter of the equivalent crystal is greater than 
the equilibrium value, the net change is still positive, but, if 
it is smaller than the equilibrium value, the net change is then 
negative. Again, this is only an approximation, and its validity 
will be tested by several applications of the method. 

Finally, the last term in the chemical energy equation, e~° s 
(see eq. (9)), is introduced to compensate for the information 
lost when breaking up the process of formation of the alloy 
from its pure constituents into two independent processes (as 
shown in fig. 6). This term is proportional to the lattice param- 
eter mismatch of the constituents. Its main role, as noted 
earlier, is to properly account for the asymptotic behavior of 
the excess energy in the formation of the alloy from its 
constituents, but it is almost of no relevance in the region of 
interest (when the typical distances in the alloy are comparable 
to the distances in the equilibrium metals). 


Determining A ab and A ba 

In order to fix the perturbative parameters A, we will use 
the experimental values for the heat of solution (in the dilute 
limit). Let's call E BA the heat of solution of an impurity B in 
a host A. This quantity is defined as 


pBA _ dAH 

d* v-o 


(19) 


where AH is the heat of formation of the compound A t _ X B X , 
x being the concentration of B atoms. Correspondingly, 


gAB _ dAH 
dx 


( 20 ) 


The heat of formation AH is defined as the difference in 
energy between the constituent atoms of the compound at 
equilibrium and their composition-averaged bulk values in the 
pure crystals. Later in this section, we will provide an 
algorithm for the calculation of AH. 

Equations (19) and (20) form a system of two coupled 
equations that are solved with respect to A ab and A ba . 
Consider a binary alloy A-B on a lattice of fixed symmetry. 
We first compute the excess internal energy AE for some 
ordered alloys in reference to its phase separation limit. For 
the ordered structures, we first consider the case of a face- 
centered-cubic (fee) lattice. If we assume that only clusters 
consisting entirely of nearest neighbors are important, then 
we only deal with structures of the form A m B 4 _ m (where the 
largest cluster of nearest neighbors is a tetrahedron as shown 
in fig. 1 1). 

These clusters are the building blocks of the corresponding 
ordered fee alloys A fn B 4 _ m (see fig. 12). The excess internal 
energy for the ordered alloy per atom is written as 


A £,„(/-) - £•„,(/•) - ^ E 4 (a?) - ^ EMc) (21) 

where a A and a B are the equilibrium lattice constants of pure 
A and B metals, and E () (a B ) and E 4 (a A ) are the cohesive 
energies E?. and E A , respectively. With these definitions, 
E„,(r) represents the binding energy curve of the ordered 
alloy. 



(fee) (L1 0 ) fee 



(L1 2 ) (L1 2 ) 

Figure 1 1.— Cluster of nearest neighbors. 
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Figure 12.— Ordered fcc, A, n B m 4 alloy structures. 

To compute A E m (r) using the method described in previous 
sections, let us consider a compound A m B A _ nr where the 
elementary structure contains N A atoms of type A and N B of 
type B, arranged in such a way that each one of the A atoms 
has N aa (M aa ) (next) nearest neighbors of type A and N AB 
(M ab ) (next) nearest neighbors of type B. The B atoms have, 
each, N B8 (M bb ) (next) nearest neighbors of type B and N BA 
(M ba ) (next) nearest neighbors of type A. The strain energy 
for each A atom is 


ef ( r) = E '\ 1 - ( 1 + a sA ) exp ( -a \ A ) ] for a\ A = 


r ~ a* 
/ 

( 22 ) 


where a A is a solution of 


NRfr exp (-a A R\) + MR?* exp - 


- 

1 

p 

tk. 

+ 

1- 



, 

V 


7 


= ^AA rl \ A exp ( ~ <*A r \ ) + N A B r \ A exp [ - (ot A + A ba ) /*| ] 


+ M AA r 2 A exp 


+ M AB n A exp 


-(«/. + — )^2 


a A + &BA + ~ ) r 2 


(25) 


for 


„ V2 <1 

R \ = —*a r 2 = a A r, = —r r 2 = r 


(26) 


where s A = + 1 if a A < 0, s A = — 1 if a A > 0, and N = 12 
(for fcc structures) and M = 6 are the number of nearest and 
next-nearest neighbors, respectively. 

A similar expression to equation (24) is found for the 
chemical energy of the B atoms just by exchanging A ~ B 
in equations (24) to (26). It is clear then that the excess energy 
(eq. (21)) is just 

A E m (r) = me A n) (r) + (4 - m)4" 0 (r) (27) 


with 


and 


el m) (r) =ef(r) + e?(r;m) 


(28) 


and for each B atom is 


elTHr) =ef(r) + e?(r;m) 


(29) 


ef(r) = Efl 1 - ( 1 + a.vfi)exp (-a f fi)] for a sB = 

(23) 

The chemical energy for each A atom is 

e?(r;m) =^£^[1 -(1 + a A ) exp(-a^)] exp(-a* /1 ) 

for a A = ^23*. (24) 


As noted before, we will determine the values of A ab and A ba 
so that our final results exactly reproduce the experimental 
heats of solution, which should be computed for the disordered 
alloys A X B j_ r . An easy way of obtaining the excess energy 
expressions for the disordered structures A X B\ _ r from those 
for the ordered ones is by following the prescription of 
Connolly and Williams (ref. 5), which takes advantage of the 
cluster expansion of Sanchez and de Fontaine (ref. 6). In this 
framework, the excess energy A E m (r) is written in terms of 
many-body interaction potentials v n (r): 

A E m (r) = ^ v„(r)i, m (30) 

n 
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TABLE I. -MULTISITE CORRELATION FUNCTIONS 
FOR FCC STRUCTURES 


where 


Formula 

Structure 

Multisite correlation function 

Co 

Cl 

^2 

Cs 

u 

A 

fee 

I 

1 

1 

1 

1 

A<B 

LU 

1 

Vi 

0 

-V4 

-1 

(AB)i 

L 1 o 

1 

0 

— */3 

0 

1 

AB, 

L l t 

1 

Vi 

0 

Vi 

-l 

B 

fee" 

l 

— 1 

1 

-1 

1 


where £ are multisite correlation functions defined on an n lh 
order cluster, given by 


^ N H °p« 

” I/M 


(31) 


where a is a spin-like variable which takes the values + 1 or 
- 1 when the lattice site p is occupied by an A or a B atom. 

The sum is over all // th -order (n = 0 4) clusters of a given 

type in the lattice, and N n is the total number of such clusters. 
For fee structures, the £'s are those shown in table I. 

Inverting equation (30), we obtain the following expression 
for the many-body potentials: 


„,,(/•)= y A£„.(')[r']„ 


(32) 


c,Jx) = Y ( 1 - 2 *)"!$ = ( m )x m ( l -X) 4 - 

n = 0 


( 35 ) 


The heat of formation of the ordered structures A E m (r m ) 
is defined as the minimum value of the excess energy A E m (r) 
(eq. (27)), which occurs for r = r,„, and the heat of formation 
of the disordered alloys AE D (r x ,x) is defined as the minimum 
value of A£ D (r,jc) (eq. (32)), which occurs for r = r v . 

For a given pair of values ( A ab , A ba ) one can evaluate the 
heat of formation of the disordered alloys for any arbitrary 
concentration x and then compare the predicted value of the 
derivatives at x = 0 and 1 with the experimental values for 
the heats of solution (See ref. 7 and eqs. (19) and (20)). This 
evaluation is, in principle, feasible but obviously impractical. 
Computing some selected heats of formation for some specific 
concentrations is clearly desirable. One choice of these selected 
concentrations is the set (x = 0, x = *4, x = '/ 2 , x = 3 4, 
x = 1). We then interpolate the intermediate values with a 
polynomial 


P{x) = ax A + bx" + cx 2 + dx - he 


(36) 


m = 0 


so that 


with 


= — 
s /?//? , , 

1 O 


14 6 4 1 

48 0-8-4 

6 0 -12 0 6 
4-8 0 8-4 

1-4 6-4 1 


(33) 


(37) 


P( 0) = / > ( 1 ) = 0 
P( 1/4) = AEp ( /'i / 4 ,* x = 1/4) = y i 
P(l/2) = A E D (r W2 ; x= 1/2) = y 2 
P( 3/4) = A E D (r 3/4 ; x = 3/4) = y 3 


Solving this system of linear equations we obtain the 
coefficients: 


The excess energy of completely disordered alloys is 
estimated by 


AE D (r,x) = Y (1 ~2x)\,(r) 

/! = 0 

4 4 

= y d- 2 *)" Y A£ m<'-)[r']„ 


/;/ = 0 


// = 0 
4 

= Y c m (x)AE„,(r) 
m=0 


(34) 


256 

a - 64v 2 ( v t + V 3 ) 

6 


b = — (9vi - 12v 2 + 7y 3 ) 


16 / 57 

c = — I — 1 3vi H y-> — 7v 3 

3 \ 4 


, 16 L 9 

a = — 3vi — y 2 + ys 
3 V ’ 4 • - 


(38) 
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Finally, the heats of solution (eqs. (19) and (20)), in terms 
of the coefficients of this polynomial, become 


TABLE III. -PERTURBATIVE PARAMETERS OBTAINED 
BY FITTING TO EXPERIMENTAL VALUES 
FOR HEATS OF SOLUTION^ 


E ba = P (0)=d 


e ab 


-/>(!)= 16 



Vi + V-* 

4 


(39) 


The parameters A ab and A fi/1 are varied until equations (39) 
are exactly satisfied. By determining the final values of the 
perturbative parameters A ab and A ba in this fashion, one 
automatically determines, by recourse to equations (30) to (35), 
the heats of formation of all the intermediate random alloys 
for arbitrary values of the concentration. This determination 
enables us to obtain the bulk properties of all the possible 
alloys. The values of the relevant ECT parameters of some 
fee metals, as well as the results obtained with our formalism 
for the possible alloys of these elements, are discussed in the 
next section. 


A 

B 

Perturbative parameters. 
A" 1 

Heats of solution, 
eV/atom 

Al B 

A.4 

e ba 

E aB 

Cu 

Ag 

-0.0321 

-0.0394 

0.39 

0.25 

Ag 

Au 

-.0311 

-.0220 

-.19 

-.16 

Cu 

Au 

- .0588 

- .05095 

-.19 

-.126 

Ni 

Cu 

.02395 

-.0131 

.09 

.03 

Al 

Cu 

-.0526 

-.0626 

-.20 

-.35 

Al 

Ag 

.0475 

-.0499 

.141 

-.166 

Al 

Au 

-.0501 

-.0853 

-.80 

-1.26 

Ag 

Pd 

-.0431 

-.02033 

-.289 

-.108 

Au 

Pd 

-.0439 

-.0348 

-.355 

-.195 

Cu 

Pd 

-.04205 

-.04795 

-.436 

-.392 

Cu 

Pt 

-.0568 

-.0444 

-.532 

-.299 

Ni 

Pd 

-.0401 

- .04665 

.057 

-.088 

Ni 

Pt 

-.0603 

-.0529 

-.282 

-.330 


a Scc ref. 7. 


Results 

We applied the formalism to several binary alloys of Ni, 
Cu, Al. Au, Ag, Pd, and Pt. As with previous calculations 
(ref. 8), in this first application we adopted the approximation 
that, below the melting point, the alloying energy is 
independent of temperature. In later efforts, we will extend 
our formalism to include temperature effects, as well as local 
relaxation effects which, for simplicity, were ignored in the 
present calculations. Local relaxation effects are, however, 
important, especially when the atomic sizes of the constituents 
differ greatly (ref. 8). 

Before applying our formalism to defect structures, local 
relaxation should be incorporated, although in the present 
report, it would just translate into a slightly different set of 
values of the perturbative parameters A ab and A Ba from the 
ones reported here. For completeness, table II shows the ECT 
parameters for the pure metals used in this calculation. 


TABLE II. -EXPERIMENTAL INPUT AND COMPUTED CONSTANTS 


Element 

P 

Screening 

length 

/ 

a, 

A 1 
(a) 

Electronic 

screening 

length. 

X 

Experimental 

cohesive 

energy, 

E 

Equilibrium 

lattice 

constant, 

a 

Al 

4 

0.336 

2.105 

0.944 

3.34 

4.04 

Cu 

6 

0.272 

2.935 

.765 

3.50 

3.61 

Ag 

8 

0.269 

3.337 

.756 

2.96 

4.08 

Au 

10 

0.236 

4.339 

.663 

3.78 

4.07 

Ni 

6 

0.270 

3.015 

.759 

4.44 

3.51 

Pd 

8 

0.237 

3.612 

.666 

3.94 

3.89 

Pt 

10 

0.237 

4.535 

.666 

5.85 

3.92 


a Sce eq. (14) 


Table III gives the values of the parameters A ab and A ba 
that one obtains when applying the formalism described in the 
previous section. As predicted by equations (39), with these 
values of A ab and A Ba , the heats of solution are exactly 
reproduced, as defined by equations (19) and (20). Note that 
although the formal definitions of equations (19) and (20) are 
exact, the estimates from experiment are not. They are obtained, 
as others have done, from a linear approximation to the slope 
at the minimum concentration available, which are often too 
high to be accurate. Our procedure allows a method of 
improving on this estimate. This will be treated in a later work. 
The excess energy curves for ordered and disordered alloys 
(A£„,(/*) and A E D (r,x), respectively) can be written as (see 
eq. (21)): 

A E m (r) = E m ( r ) + A* 

A E D (r,x) = E x (r) -P A* (40) 

with 

+ (i -f)E? <41, 

where £,„(/*) ( £ v (r )) is the binding energy curve for the 
ordered (disordered) compound A tn B 4 _ m (A X B\^ X ). These 
functions can be accurately described with analytical forms as 

E m (r) = -£"'( 1 + aje< 

E x ( r ) = -E?( 1 +a x )e~ a 'x (42) 


n 


TABLE IV.-ECT RESULTS FOR ORDERED AND DISORDERED ALLOYS 


Formula 

Ordered alloys, A m B A _, n 

Disordered alloys, A X B\_ V 

Experiment values 

Predicted 

ordering 

energy, 

E ord- 

eV 

Lattice 

parameter, 

r nr 

A 

Heat 

formation, 

(*///)' 

eV 

Energy 

reference, 

A®, 

eV 

Cohesive 

energy. 

eV 

Lattice 

parameter, 

r v , 

A 

Heat of 
formation. 
AE D (r x ,x), 
eV 

Energy 

reference, 

A®. 

eV 

Cohesive 
energy, 
E* (exp), 
eV 

Heat of 
formation, 

AW cxp . 

eV 

Cohesive 

energy. 

eV 

CuAu 3 

4.003 

-0.0588 

3.710 

3.769 

4.008 

-0.0217 

3.710 

3.732 

-0.0311 

3.741 

0.0371 

(CuAub 

3.909 

-.0818 

3.640 

3.722 

3.918 

-.0308 

3.640 

3.671 

- .0539 

3.693 

0.051 

Cu 3 Au 

3.771 

-.1646 

3.570 

3.735 

3.798 

-.0211 

3.570 

3.598 

-.0452 

3.615 

.1369 

CuAg 3 

3.992 

0.0366 

3.095 

3.058 

3.996 

0.0537 

3.095 

3.041 



0.0171 

(CuAgb 

3.885 

.0565 

3.230 

3.173 

3.890 

0.0803 

3.230 

3.150 

— 

— 

.0238 

Cu 2 Ag 

3.756 

.0187 

3.365 

3.346 

3.765 

.0668 

3.365 

3.298 

— 

— 

.0488 

AICu 3 

3.688 

-0.1083 

3.460 

3.568 

3.698 

-0.0520 

3.460 

3.512 

-0.0728 

3.533 

0.0563 

(AICu ) 2 

3.788 

-.1012 

3.420 

3.521 

3.799 

- .0569 

3.420 

3.477 

- .0938 

3.514 

.0443 

AI 3 Cu 

3.897 

- .0999 

3.380 

3.480 

3.915 

-.0378 

3.380 

3.418 

-.0521 

3.432 

.0621 

NiCu 3 

3.590 

0.0052 

3.734 

3.729 

3.592 

0.0083 

3.734 

3.726 

0.0093 

3.725 

0.0031 

(NiCu) : 

3.568 

.0171 

3.967 

3.950 

3.568 

.0147 

3.967 

3.952 

.0184 

3.949 

-.0024 

Ni 3 Cu 

3.548 

.0200 

4.201 

4.181 

3.548 

.0140 

4.201 

4.187 

.0179 

4.183 

-.006 

N i Ft 3 

3.857 

-0.1129 

5.496 

5.609 

3.863 

-0.0587 

5.496 

5.555 

-0.0641 

5.560 

0.0542 

(NiPt ) 2 

3.782 

-.1390 

5.142 

5.281 

3.790 

-.0752 

5.142 

5.217 

-.0960 

5.238 

.0638 

N i 3 Pt " 

3.669 

-.1958 

4.789 

4.985 

3.689 

-.0545 

4.789 

4.843 

- .0695 

4.858 

.1413 

NiPd 3 

3.817 

-0.0493 

4.064 

4.113 

3.822 

-0.0088 

4.064 

4.073 

-0.0123 

4.076 

0.0405 

(NiPd ) 2 

3.738 

-.0411 

4.187 

4.228 

3.744 

-.0022 

4.187 

4.189 

-.0055 

4.192 

.0389 

Ni 3 Pd 

3.639 

-.0581 

4.311 

4.369 

3.649 

.0048 

4.311 

4.306 

.0050 

4.306 

.0629 

CuPt 3 

3.877 

-0.0891 

5.262 

5.351 

3.879 

-0.0641 

5.262 

5.326 

-0.0723 

5.334 

0.025 

(CuPt ) 2 

3.818 

-.1502 

4.675 

4.825 

3.825 

-.0986 

4.675 

4.774 

-.1149 

4.790 

.0516 

Cu 3 Pt " 

3.729 

-.1970 

4.087 

4.284 

3.747 

-.0859 

4.087 

4.173 

-.0977 

4.185 

.1111 

CuPd 3 

3.830 

-0.1040 

3.830 

3.934 

3.837 

-0.0687 

3.830 

3.899 

-0.0818 

3.912 

0.0453 

(CuPdb 

3.771 

-.1306 

3.720 

3.851 

3.777 

-.0912 

3.720 

3.811 

-.1109 

3.831 

.0394 

Cu 3 Pd " 

3.693 

-.1387 

3.610 

3.749 

3.706 

-.0728 

3.610 

3.683 

- .0926 

3.703 

.0659 

AuPd 3 

3.933 

-0.0597 

3.900 

3.960 

3.935 

-0.0429 

3.900 

3.943 

-0.0476 

3.948 

0.0168 

(AuPd ) 2 

3.975 

-.0966 

3.860 

3.957 

3.979 

- .0667 

3.860 

3.927 

- .0807 

3.941 

.0299 

Au 3 Pd ’ 

4.018 

-.1061 

3.820 

3.926 

4.026 

-.0579 

3.820 

3.878 

-.0746 

3.895 

.0482 

AgPd 3 

3.921 

-0.0330 

3.695 

3.728 

3.921 

-0.0303 

3.695 

3.725 

-0.0273 

3.722 

0.0030 

(AgPd) : 

3.954 

-.0761 

3.450 

3.526 

3.959 

-.0523 

3.450 

3.502 

-.0520 

3.502 

.0238 

Ag 3 Pd 

4.000 

- .0956 

3.205 

3.301 

4.011 

-.0471 

3.205 

3.252 

-.0524 

3.257 

.0485 

AIAu 3 

4.075 

-0.3155 

3.463 

3.778 

4.075 

-0.2068 

3.463 

3.670 

-0.2723 

3.735 

0.1087 

(AIAu ) 2 

4.070 

- .3054 

3.316 

3.621 

4.070 

-.2434 

3.316 

3.559 

-.3499 

3.666 

.0620 

AI 3 Au 

4.062 

-.2006 

3.286 

3.487 

4.062 

-.1636 

3.286 

3.450 

-.2095 

3.495 

.037 

AlAg, 

4.072 

-0.0420 

3.055 

3.097 

4.077 

-0.0160 

3.055 

3.071 




0.0260 

(AIAg) : 

4.074 

.(XXX) 

3.15 

3.15 

4.073 

-.0018 

3.150 

3.152 



-.0018 

Al,Ag 

4.073 

.0342 

3.245 

3.21 1 

4.067 

.0127 

3.245 

3.232 

— 

— 

-.0215 

AgAu 

4.079 

-0.0401 

3.535 

3.575 

4.079 

-0.0307 

3.535 

3.566 

-0.0341 

3.569 

0.0094 

(AgAu) 

4.081 

- .0550 

3.315 

3.370 

4.081 

-.0425 

3.315 

3.357 

-.0482 

3.363 

.0125 

Ag 3 Au 

4.083 

- .0475 

3.117 

3.165 

4.083 

-.0335 

3.117 

3.150 

- .0372 

3.154 

.0140 


| 
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where 

r - r m . r — r x 

a m = and a x = 

/ 

l m f .i- 

In the equations E"' (E c x ) is thus the cohesive energy of the 
ordered (disordered) alloy, and r m (r x ) is the equilibrium 
lattice parameter. This result is in agreement with previous 
attempts to represent the binding energy curves for alloys. 
Moreover, equation (34) implies that 

E x (r) = £ c m U)E m {r) (43) 

m 

and 

A.r= X) c '» (v)A ” (44) 

m 

so that only the knowledge of E”\ a nr /,„, and A" is required, 
as the quantities E x , r x , l x , and A“ can be obtained from 
equations (43) and (44). These results are listed in table IV 
for some of the binary alloys studied. Table IV also gives the 
predicted values for the heats of formation of ordered and 
disordered structures, as well as the experimental values. 


Finally, table IV includes the ordering energy, defined as 
EW = AE d (r x .x) - AE„,(r m ) (45) 

for x = ml 4 and m = 0, 1,2,3, 4. 

In figure 13 we show a particular case, the Ll 0 structure 
of the Ag-Pd alloy. This figure highlights the relative 
contributions of the strain and chemical energies to the heat 
of formation A E 2 (r). Figure 14 compares the experimental 
values of the heat of formation (when available) with the 
predicted values obtained with our formalism. In most cases, 
especially for those alloys where the lattice mismatch is not 
large, the agreement is excellent. Good results are obtained 
in general for all the cases studied, which include two liquids: 
Al-Cu and Al-Au. Figure 15 compares the experimental values 
of the cohesive energies of the disordered compounds with 
the predicted values, both in reference to the values one would 
obtain if Vegard’s law (linear average of the pure metals 
cohesive energies) were valid. In all cases the agreement is 
excellent, even in those situations where there is noticeable 
departure from Vegard’s law. As a final test, we compare the 
bulk modulus predicted by ECT with experiment and the first- 
principles calculations of Wei et al. (ref 9) for the ordered 
compounds of Cu and Au. The results are given in table V. 
This is a severe test of the predictions since the bulk modulus 
is related to the second derivative of the binding energy curve. 
Again, the agreement is excellent and comparable to the first- 
principles results. 



Figure 13.— Excess energy A E{r) of ordered compound (AgPd) 2 as function of lattice parameter r also showing the two contributions to excess energy 

Fslr;iin!*^ F c h cm (r). 
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Figure 14.— Heat of formation as function of concentration 
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Figure 14.— Concluded. 
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Figure 15.— Cohesive energy as function of concentration for various alloys. 
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Cohesive energy, E 





Figure 15.— Concluded. 
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TABLE V. -COMPARISON OF EXPERIMENTAL RESULTS. ECT 
PREDICTIONS. AND FIRST-PRINCIPLES 
CALCULATIONS (FPPL) 


Formula 

Type of result 

Lattice 

parameter. 

Cl. 

A 

Screening 

length, 

/. 

A 

Cohesive 

energy, 

£,. 

eV 

Bulk 

modulus, 
B . 

GPa 

Cu 

Experiment 

3.615 

0.276 

3.49 

138 


ECT predictions 

3.615 

.272 

3.50 

142 


FPPL 

3.577 

.302 

4.33 

144 

Cu^Au 

Experiment 

3.743 

0.267 

3.64 

148 


ECT predictions 

3.771 

.255 

3.735 

165 


FPPL 

3.738 

.301 

4.37 

140 

(CuAu) 2 

Experiment 

3.876 

0.254 

3.74 

163 


ECT predictions 

3.909 

.246 

3.722 

171 


FPPL 

3.887 

.275 

4.40 

162 

CuAu^ 

Experiment 

3.982 

0.247 

3.79 

170 


ECT predictions 

4.003 

.240 

3.769 

177 


FPPL 

3.991 

.248 

4.37 

194 

Au 

Experiment 

4.078 

0.244 

3.81 

171 


ECT predictions 

4.078 

.236 

3.78 

181 


FPPL 

4.106 

.253 

4.35 

180 


Discussion 

The results described in the previous section show how our 
calculations compare with the body of experimental data 
available, from which we extracted the information (i.e., heats 
of solution) needed to determine the parameters included in 
our theory. 

In what follows, we will show how our results compare with 
some recent theoretical studies which include first-principles 
total energy calculations, based on solving the three-dimensional, 
nearest-neighbor, fee lsing model with volume-dependent 
interaction energies, local density approximation (LDA) calcu- 
lations, and embedded atom method (EAM) results (refs. 9 
to 11). 

Comparison with First-Principles Calculations 

Figure 16 shows the first of these comparisons for the 
particular case of Cu-Au alloys. All three sets of results for 
the binding energy curves of the ordered compounds A n fi A ~ m 
(i.e., experimental values (EXP, ref. 7), first-principles 
(FPPL, ref. 9), and ECT) can be represented by simple 
Rydberg functions of the form 

E m (r) = -£;”( 1 +a m )e-** (46) 

with 

a,„ = - (r ~ rj 

hn 



Lattice parameter, A 

Figure 16.— Binding energy as function of lattice parameter for gold (Au). 
(See ref. 9.) 


where q = (3/167t) 1/3 for fee only. Table V gives the values 
of E"\ r m , and l m for all three sets of results (EXP, FPPL, 
and ECT). For clarity, we first show the binding energy curves 
for gold in figure 16 as obtained with the parameters of 
table V. Figure 17 expands on these results showing the 
corresponding results for Cu, Au, and the three intermediate 
ordered compounds. 
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Figure 17.— Binding energy as function of lattice parameter for ordered 
compounds Cu m Au 4 _ m . 


Comparison with Other Semiempirical Methods. 

In a recent study (ref. 11), the embedded atom method 
(EAM) was applied to the calculation of formation enthalpies 
and lattice parameters of Pd-Ni alloys. The embedding func- 
tions used were those determined for the pure metals, and the 
original EAM prescription was used. The comparison with 
experimental results is disappointing. Since then, EAM has 
been reformulated for alloys in order to improve the accuracy 
of its predictions. 

■ Experiment 
ECT results 



Figure 1 8. — Heat of formation as function of concentration for Pd-Ni alloys. 



Figure 19.— Deviation from VegarcTs law for lattice parameters of Pd-Ni alloys 

as function of concentration. 

For comparison, figure 18 shows the ECT values, the EAM 
results, the experimental values, and the predictions of 
Miedema's empirical model for the formation enthalpy of the 
Pd-Ni alloys as a function of composition. Figure 19 shows 
the corresponding results for the lattice parameters. 

Decoupling the Strain and the Chemical Energy 

As pointed out before, our formalism is based on the 
assumption that the structural changes in alloy formation (strain 
energy) and the composition changes (chemical energy) can 
be treated separately. This is later corrected (in an approximate 
way) by introducing a geometrical factor in the expression for 
the chemical energy (see eq. (9)). It was later argued that in 
spite of the fact that this term is strictly necessary in order 
to ensure the correct asymptotic behavior of the formation 
enthalpies, it is not fundamentally important for the study of 
the alloy properties near equilibrium. 

Not including this term in the chemical energy would 
effectively decouple the chemical and the strain energies 
leading to an approximation to what is already an approximate 
method. However, we will show below that this approximation 
leads to an interesting result concerning general properties of 
alloys, which in most cases can make a full calculation of alloy 
structure with any available method, unnecessary. 

If the geometry factor in the expression for the chemical 
energy is left out (e~° s in eq. (9)), the chemical energy is 
dependent only on chemical composition and the strain energy 
carries all the information about the atomic distribution and 
composition in the alloy. The values of the parameters A ab 
and A ba obtained are then slightly different from the ones 
quoted earlier, as shown in table VI, where the two sets of 
A’s are listed. The ensuing results for the enthalpies of 
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TABLE VI. -COMPARISON OF COUPLED 
AND UNCOUPLED VALUES OF THE 
PERTURBATIVE PARAMETERS 


A 

A 

A* B 

^ BA 

Uncoupled 

Coupled 

Uncoupled 

Coupled 

Cu 

Ag 

-0.0380 

-0.0321 

-0.0289 

-0.0394 

Ag 

Au 

-.0311 

-.0311 

- .0220 

-.0220 

Cu 

Au 

-.0711 

-.0588 

-.0404 

- .05095 

Ni 

Cu 

.0260 

.02395 

-.0130 

-.0131 

A1 

Cu 

-.0415 

-.0526 

-.0775 

-.0626 

Al 

Ag 

.0487 

.0475 

- .0492 

-.0499 

Al 

Au 

-.0501 

-.0501 

- .0853 

-.0863 

Ag 

Pd 

-.0381 

-.0431 

-.0246 

-.02033 

Au 

Pd 

-.0385 

-.0439 

-.0411 

-.0348 

Cu 

Pd 

- .0488 

-.04205 

-.0413 

-.04795 

Cu 

Pt 

-.0653 

-.0568 

-.0367 

-.0444 

•Ni 

Pd 

-.0491 

-.0401 

-.0376 

-.04665 

Ni 

Pt 

-.0719 

-.0603 

-.0431 

-.0529 



formation are still good and in some cases indistinguishable 
from the previous results. Figure 20 displays, for comparison, 
the original results, the new, approximate results, and the 
corresponding experimental values for some of the compounds 
studied. 

As before, the results obtained for the binding energy curves 
of the disordered alloys, as well as the corresponding ones 
for the ordered alloys, can be accurately represented by simple 
analytical expressions: 

A £,„(/•) = - £<"'( 1 + + A”; a' m - ~ (/- - r,„) 

*/« 

(47) 

A£ 0 (r,x) = - E'( 1 + aj)e“"« + A“- a' v = j (r - r v ) 

.V 

(48) 




Figure 20.— Comparison of approximate ECT results and results shown in figure 14. 


j 



Ni concentration 




Figure 20.— Concluded. 


where the quantities E"\ r m , /,„, and A“, as well as E x c , r x , 
/,, and A,°°, are also different from the ones one obtains when 
applying the original formalism. However, in this case, the 
sets of values r m , /„,) and (£ ( v , r x , l x ) are the same when 
x = ml 4. This means that, apart from the reference energies 
A“ and A", the binding energy curves for the disordered and 
ordered alloys of the same concentration coincide exactly. The 
only difference arises in the reference energy term. 

If we insert equations (47) and (48) into equation (34), we 
find that two separate conditions are then satisfied: 

4 

A.r= Y. c m (x)*n, (49) 

m ~ 0 


and 

4 

£,'(1 + a' x ) exp(-a“) = ^ c m (x)E"'( 1 + a m ) exp 

m=0 

(50) 

These results were also valid before, but the fact that the 
sets of values (E"\ r m , l m ) and (E'. t r x , l x ) are identical when 
x = ml 4 we may write equation (50) in the following fashion, 
reflecting the equivalence of the binding energy curves for 
specific concentrations (i.e., E x (r) = E,„(r) for x = m/4): 
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E„(r) = Yj c >« ( x ~f) E m(r); n= 0 4 (51) 

m = 0 ' 4 ' 

Equation (51) is just a linear system of equations for £,(/* ), 
which yields the result 

( m\ m 

1 - E a (r) +-E 4 (r) (52) 

This result, which is only valid within the context of this 
approximation (when the strain and chemical energy are 
“decoupled”) states that Vegard’s law is valid for the binding 
energy curve and not just for the individual parameters in them. 
By virtue of equation (34), we can extend this result to the 
disordered alloys, obtaining 

E D (r,x) = (1 -x)E 0 (r) + xE 4 (r) (53) 

With these results, certain alloy properties can then be obtained 
without a full many-atom calculation. One would expect these 
predictions to be valid according to the relevance of the 
coupling between the chemical and strain energies. Figure 20 
shows that, in some cases, this is an excellent approximation. 
The prescription for its use, therefore, is to use the UBER 
to generate the binding energy curves for the pure materials 
and then to build the alloy curve as prescribed above. 

Note that this approximation is limited, not only by the 
underlying assumptions, but also by its inability (by definition) 
to reproduce the correct asymptotic behavior of the excess 
energy. The range of validity is thus restricted to the “local” 
regime, that is, when the distances involved are comparable 
to those in the pure metals near equilibrium. 

Comparison with Other First-Principles Calculation— Test 
of the Sum Rule for Binding Energy Curves 

In a recent calculation Terakura et al. (ref. 10) analyzed the 
phase stability of several binary alloy systems composed of 
Ni, Pd, or Pt as one element and Cu, Ag, or Au as the other. 
From their LDA band calculations they found that their pre- 
dictions for the excess energy A E,„(r) admitted a simple 
parameterization of the form 



where the parameters p m , q m and r m are those listed in table 
VII. They also found that a good fit was obtained for n = 3.5. 

The objective in this section is to show that even first- 
principles results, to a certain extent, follow the general sum 
rule derived in the previous section (eq. (52)). To that effect. 


TABLE VII. — FITTING 
PARAMETERS OF EQUATION (45) 

[Ref 10.] 


System 

P 

<7 

r 

Au 

8.584 

11.677 

2.155 

CuAu 3 

8.217 

11.013 

1.932 

(CuAu) 2 

7.877 

10.450 

1.788 

Cu 3 Au 

7.474 

9.749 

1.586 

Cu 

7.055 

9.033 

1.410 

Ag 

8.075 

10.444 

1.513 

CuAg, 

7.819 

10.076 

1.502 

(CuAg) 2 

7.588 

9.787 

1.519 

Cu } Ag 

7.322 

9.406 

1.469 

AgAu 3 

8.447 

11.347 

1.960 

(AgAu) 2 

8.331 

1 1 .073 

1.815 

Ag 3 Au 

8.207 

10.769 

1.662 


we note that in order for this rule to hold true, equation (51), 
together with equations (49) and (50), must be satisfied for 
any arbitrary choice of the parameterization of the binding 
energy dependence on the lattice parameter. Introducing 
equation (54) into equation (34), we obtain the following 
expression for the excess energy of the corresponding 
disordered compounds: 

AE D (r,x) = % - =7 + A, (55) 

r r 

where 

4 

p x" = Yi c <» (x) p»" 

w = 0 
4 

q::= y (56) 

w= 0 
4 

Cm(x)r , m 

m = 0 

If the first-principles results behave according to equation 
(52), then a direct consequence of the sum rule (eq. (34)), 
namely, 

m\ 

r t x = — j — AE,,,^) m = 0,...,4 (57) 

would also have to be satisfied. In terms of the coefficients 
used in Terakura’s work, this condition can be written as a 
set of simultaneous requirements on these coefficients: 
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TABLE VIII -COMPARISON OF THE SETS 
OF PARAMETERS (P.Q.R) AND (p,q,r) 
AS DEFINED BY EQUATION (57) 


System 

P/p 

Q/q 

R/r 

AgAu 3 

(AgAu) 2 

Ag 3 Au 

8.46/8.45 

8.33/8.33 

8.21/8.21 

11.37/11.35 

11.07/11.07 

10.77/10.77 

1.98/1.96 

1.82/1.82 

1.66/1.66 

CuAg, 

(CuAg), 

Cu,Ag 

7.84/7.82 

7.60/7.59 

7.34/7.32 

10.11/10.08 

9.77/9.79 

9.42/9.41 

1.51/1.50 

1.50/1.52 

1.46/1.47 

CuAu 3 

(CuAu) 2 

Cu 3 Au 

8.26/8.22 

7.91/7.88 

7.52/7.47 

11.07/11.01 

10.45/10.45 

9.78/9.75 

1.95/1.93 
I. 77/1. 79 
1.59/1.59 


P 2 % = P\ 

II 


oi 

^ 1/4 — r 1 

P 2 U2=P2 

Q'l/2 — #2 

R\/ 2 — 

Pl% = Pi 

Q"/ 4 = <?3 

^3/4 = r 3 


The results are shown in table VIII. Although the comparison 
is not exact (i.e., the (P,Q,R) and (p,q,r) coefficients are not 
exactly the same), the agreement is surprisingly good, making 
the use of equation (52) a viable alternative to Vegards’s law 
for obtaining a more accurate and complete calculation of 
general alloy properties. 

Conclusion 

The equivalent crystal theory, originally developed for the 
study of pure metals and semiconductors, has been extended 
to include binary alloys. A simple formalism, inspired by the 
concept of equivalent crystals, was derived and, with minimum 
experimental input, applied to a large number of metallic alloys 
of Cu, Ni, Pd, Pt, Au, Al, and Ag. The results, which in all 
cases closely follow the experimental measurements of certain 


properties of the alloys for all ranges of compositions, compare 
favorable with other approaches. By construction, the method 
allows for a simple treatment of defects in alloys as well as 
extensions to include temperature effects , making this new 
equivalent crystal theory of alloys a versatile tool for several 
applications of interest. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 4, 1991 
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